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EXECUTIVE  SUMMARY 


OBJECTIVE 

Predict  the  polarization  dependence  of  infrared  sea  radiance  (brightness). 

APPROACH 

The  Fresnel  coefficients  are  applied  to  the  reflection  of  unpolarized  monochromatic  radiation  by  a 
sea  water  capillary  wave  facet.  The  result  for  a  single  tilted  facet  is  integrated  over  all  possible  tilts 
and  wave  numbers  to  predict  the  polarization  dependence  of  radiance  received  from  the  wind-ruffled 
ocean  surface. 

RESULTS 

Equations  valid  in  the  infrared  are  derived  for  horizontal,  vertical,  and  unpolarized  sea  radiance. 
They  show  the  following  regularities:  (1)  horizontal  and  vertical  sea  reflectivity  sum  to  the  unpolar¬ 
ized  sea  reflectivity,  (2)  horizontal  and  vertical  sea  emissivity  sum  to  the  unpolarized  sea  emissivity, 
and  (3)  horizontal  and  vertical  sea  radiance  sum  to  the  unpolanzed  sea  radiance. 
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1.  The  coordinate  system.  The  tilted  facet  passes  through  the  origin  at  the  point  of 
reflection.  The  X-Y plane  is  horizontal.  Unit  vector  Ug  points  to  the  source. 

Unit  vector  U,-,  coincident  with  the  optical  axis  of  the  polarizer,  points  to  the 
receiver  . 

2.  Incident  coherent  amplitudes.  The  view  is  along  vector  Ug  toward  the  source 

3.  Reflected  coherent  amplitudes.  The  view  is  against  vector  U,-  toward  the  facet. 
All  vectors  shown  here  are  normal  to  U;-  in  the  plane  of  the  paper 

4.  Sea  water  Fresnel  coefficients  at  a  wave  number  of  1000  cm~‘' 

(a  wavelength  of  10  jim) . 

5.  Sea  water  Fresnel  coefficients  at  a  wave  number  of  2500  cm"^ 

(a  wavelength  of  4  pm) . 
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6.  Sea  radiance  in  the  long-wave  band  for  a  rough  sea.  The  wind  speed  is  10  m  s■^ 
the  solar  zenith  angle  is  70  degrees,  and  the  view  is  directly  up  wind  from  a  height 
of  20  m  looking  down  at  an  angle  of  1  degree  with  respect  to  the  horizon.  Numbers 

in  parentheses  refer  to  the  polarization  of  each  component .  17 

7.  Sea  radiance  in  the  mid-wave  band  for  a  rough  sea.  The  geometry  is  the  same 

as  for  the  previous  figure.  Numbers  in  parentheses  refer  to  the  polarization  of  each 
component .  18 

8.  Sea  polarization  in  the  mid-  and  long-wave  bands  for  a  rough  sea.  The 

conditions  are  the  same  as  for  figure  6.  A  negative  degree  of  polarization 
means  that  emission  dominates  reflection;  a  positive  degree  of  polarization 
means  that  reflection  dominates  emission.  Numbers  in  parentheses  refer  to  the 
equivalent  black  body  temperature  of  the  corresponding  signal  .  19 
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INTRODUCTION 


When  a  linearly  polarized  plane  wave  of  electromagnetic  radiation  is  incident  on  the  boundary 
between  two  homogeneous  materials  with  different  optical  properties,  it  is  well  known  that  the 
reflected  plane  wave  is  described  by  the  Fresnel  equations  (Bom  &  Wolf,  1970,  and  Stratton,  1941). 
The  Fresnel  equations  arise  from  requiring  a  match  between  the  components  of  Maxwell’s  equations 
on  either  side  of  the  boundary.  The  Fresnel  equations  show  that  the  amplitudes  and  phases  of  the 
reflected  waves  are  different  for  waves  incident  with  electric  vectors  parallel  (“p”)  and  perpendicular 
(“s”)  to  the  plane  of  incidence  ^  and,  furthermore,  that  these  two  components  are  mathematically  sep¬ 
arable  and  therefore  independent  of  one  another.  When  the  two  components  are  recombined,  differ¬ 
ences  that  have  arisen  between  them  during  reflection  lead  to  elliptically  polarized  light.  Electromag¬ 
netic  theory  thus  shows  that  a  surface  introduces  a  change  in  the  state  of  polarization  of  electro¬ 
magnetic  radiation. 


THE  PROBLEM 


The  purpose  of  this  report  is  to  explain  how  the  equations  for  Fresnel  reflection  can  be  used  to  pre¬ 
dict  the  polarization  of  sea  radiance.  The  experiment  we  wish  to  consider  is  as  follows;  (1)  Observe 
the  radiance  (brightness)  of  the  sea  with  an  instrument  that  responds  in  the  infrared  region  of  the 
electromagnetic  spectrum.  (2)  Provide  for  the  insertion  and  removal  of  a  polarizer  in  front  of  the 
radiometer.  If  the  polarizer  is  perfect,  it  will  completely  transmit  radiation  whose  electric  vector  is 
parallel  to  its  axis  of  polarization,  and  the  polarizer  will  completely  block  radiation  whose  electric 
vector  is  perpendicular  to  its  axis  of  polarization.  (3)  Point  the  instrument  toward  the  sea  and  make 
three  measurements:  A),,  the  radiance  with  the  polarizer  axis  horizontal;  N^,  the  radiance  with  the 
polarizer  axis  vertical;  and  N,  the  radiance  with  the  polarizer  removed. 

The  problem  is  to  predict  the  values  of  Nh,  N^,  and  N.  The  solution  to  this  problem  has  already 
been  given  (Gregoris  et  al.,  1992),  and  comparisons  have  been  made  with  data  (Cooper  et  al.,  1992). 
In  this  report,  we  supply  an  analysis  to  clarify  the  physics  involved. 


’  The  plane  of  incidence  contains  the  incident  ray,  the  boundary  normal,  and  the  reflected  ray. 
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THE  COORDINATE  SYSTEM 


We  will  analyze  the  problem  in  successive  stages  using  the  notation  of  Stratton  (1941)  as  much  as 
possible.  First,  we  consider  the  reflection  of  monochromatic  radiation  from  a  flat  sea  water  surface 
tilted  with  respect  to  the  horizon.  This  case  represents  the  situation  that  is  assumed  to  occur  when  a 
single  capillary  wave  facet  whose  radius  of  curvature  is  much  greater  than  the  wavelength  of  the 
radiation  is  exposed  to  the  sky  or  sun  (Walker,  1994  and  Zeisse,  1995a  &  1995b).  Figure  1  shows  the 
geometry  of  reflection.  Since  the  distribution  of  capillary  wave  slopes  depends  somewhat  on  the 
wind  direction  it  will  prove  to  be  convenient  for  later  analysis  if  the  coordinate  system  has  a  definite 
orientation  with  respect  to  the  azimuth  of  the  wind.  A  coordinate  system  has  been  chosen  whose  ori¬ 
gin  is  the  point  of  reflection  with  the  X-axis  pointing  up  wind,  the  T-axis  pointing  cross  wind,  and  the 
Z-axis  pointing  toward  the  zenith,  such  that  a  right-handed  system  is  formed.  The  X-T  plane  is  hori¬ 
zontal  at  the  point  of  reflection.  The  tilted  facet  passes  through  the  origin. 


z 


Figure  1 .  The  coordinate  system.  The  tilted  facet  passes  through  the  origin  at  the  point  of 
reflection.  The  X-Vplane  is  horizontal.  Unit  vector  Ug  points  to  the  source.  Unit  vector 
U;.,  coincident  with  the  optical  axis  of  the  polarizer,  points  to  the  receiver. 

Define  a  unit  vector  U  =  (0,  ([))  with  polar  coordinates  6  as  the  zenith  angle,  and  (|)  as  the  azimuth. 

If  we  denote  the  Cartesian  coordinates  of  U  by  (a,  b,  c),  then  we  have 

a  =  sin  0  cos  (|) 

Z?  =  sin  0  sin  (j)  (1) 

c  =  cos0 
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Two  such  unit  vectors  are  shown  in  figure  1;  U^,  pointing  from  the  origin  to  the  source,  and  U^-, 
pointing  from  the  origin  to  the  receiver.  A  third  unit  vector,  is  normal  to  the  facet  at  the  point  of 
reflection  but  was  removed  from  the  figure  for  clarity^.  When  there  is  a  specular  reflection  at  the 
facet,  then 


Uj +U,.  =2cos  (oU„  (2) 

where  co  is  the  angle  of  incidence  and  the  angle  of  reflection. 

The  facet  slope  in  the  upwind  direction,  is  given  by  the  slope  of  the  line  formed  at  the  intersec¬ 
tion  of  the  facet  with  the  X-Z  plane.  The  facet  slope  in  the  cross-wind  direction,  is  given  by  the 

slope  of  the  line  formed  at  the  intersection  of  the  facet  with  the  Y-Z  plane.  In  Cartesian  coordinates 
these  slopes  are 


(3) 


^The  zenith  angle  of  U„  is  the  same  as  the  tilt  of  the  facet.  The  tilt  is  the  angle  of  the  steepest  ascent  within  the  facet. 
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INCIDENT  AMPLITUDES 


We  will  imagine  initially  that  the  incident  electromagnetic  wave  arises  from  a  single  radiating 
atom  in  the  source  so  that  the  incident  wave  front  is  linearly  polarized  and  perfectly  coherent.  The 
incident  E  vector,  E^,  will  have  a  unique  direction  in  this  instance,  and  we  denote  by  a  the  angle  of 
arrival,  the  angle  that  E;,  makes  with  the  plane  of  incidence  (see  figure  2).  The  incident  vector  can  be 
resolved  into  two  components:  perpendicular  to  the  plane  of  incidence  and  parallel  to  the 

plane  of  incidence.  These  incident  amplitudes  are  shown  in  figure  2  where  it  is  apparent  that 


^«  =  ^.sina 
E,p  =  E^  cos  a 


(4) 


Figure  2.  Incident  coherent  amplitudes.  The  view 
is  along  vector  Ug  toward  the  source. 
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REFLECTED  AMPLITUDES 

Next,  the  Fresnel  equations  are  applied  to  each  incident  component  to  give  the  reflected  ampli¬ 
tudes  and  phases  (figure  3): 


£’2.  =  P.exp(-/5,)£^^ 

^2;,  =  P,exp(-f6p)£„^ 

We  assume  that  sea  water  is  a  conducting  medium  for  which  6^  hp.  Equations  for  p,,  6^,  pp,  and  hp, 
are  given  in  terms  of  the  optical  constants  of  air  (medium  2)  and  water  (medium  1)  by  equations 
derived  in  Stratton  (1941).  These  equations  are  reproduced  as  (A-9)  through  (A-11)  in  Appendix  A. 


PLANE  OF 
INCIDENCE 


Figure  3.  Reflected  coherent  amplitudes.  The  view  is  against  vector 
Ur  toward  the  facet.  All  vectors  shown  here  are  normal  to  in  the 
plane  of  the  paper. 

The  polarizer  has  an  optical  axis  along  which  the  radiation  passes  and,  normal  to  this  optical  axis, 
a  polarization  axis  that  is  the  axis  of  easy  passage”  for  the  E  vector.  Let  the  optical  axis  of  the  polar¬ 
izer  coincide  with  the  direction  of  in  figure  1  and  set  the  polarization  axis  either  normal  to  the 
plane  of  vision^  (along  the  unit  vector  U/,  in  figure  3)  or  parallel  to  the  plane  of  vision  (along  the  unit 
vector  Uv  in  figure  3).  In  each  position  only  those  components  of  £2^  and  E2p  parallel  to  the  polariza¬ 
tion  axis  are  transmitted  along  the  optical  axis.  The  horizontal,  vertical,  and  unpolarized  amplitudes 
will  be  given  respectively  by 


^  The  plane  of  vision  (Gregoris,  1992)  contains  the  vectors  and  k,  where  k  is  a  unit  vector  along  the  Z  axis. 
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(6) 


^ih  =  +  ^2.cosp  +  £2p  sinp 
^2v  =-^2.  sinp  +  fjpCosp 

^2  =  +  ^2, +E2p 

In  Appendix  C,  it  is  shown  that  p,  the  angle  between  the  plane  of  incidence  and  the  plane  of  vision, 
is  given  by 


sin  p  = 


sin(<l),-(l)Jsine„ 

sino) 


(7) 


where  <^?r  is  the  azimuth  of  the  receiver  ray  (figure  1),  ^  is  the  azimuth  of  the  facet  normal,  0„  is  the 
zenith  angle  of  the  facet  normal,  and  co  is  the  angle  of  incidence. 
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COHERENT  REFLECTED  INTENSITIES 


Now  that  we  have  expressions  in  (6)  for  the  reflected  amplitudes,  the  intensities  are  found  by  com¬ 
puting  the  squared  modulus  £*■£■  for  each  case.  Applying  this  to  the  amplitudes  given  in  (6)  and 
using  (5)  and  (4)  we  find  the  corresponding  intensities 

=[i?Xsin  acospf  4-(l/2)p^pp  sin2a  sin2Pcos  6  -i-/?p(cos  a  sin  pf  Je* 

K -K  =[^.(sin  asinpf-(l/2)p^p^  sin2a  sin2Pcos  b+R^{cos  a  cos  P/]eJ  (8) 

K-Ei  =[^j(sin  af  +  p,pp  sin2acos  6+E^(cosaf]£j  • 

In  equations  (8),  =  pp2  and  6  =  ^  -  5^  are  given  by  equations  (A-9)  through  (A-12).  Rj 

and  Rp  are  the  Fresnel  reflection  coefficients.  They  are  shown  for  sea  water  in  figures  4  and  5  for 
wave  numbers  of  1000  cm“^  (10  pm)  and  2500  cm“*  (4  fim),  respectively.  Appendix  B  contains  the 
FORTRAN  program  used  for  the  calculation  shown  in  figures  4  and  5.  The  program,  taken  from  the 
equations  of  Stratton  and  employing  his  notation  for  the  most  part,  is  a  simplified  version  of  the 
complete  treatment  (summarized  in  Appendix  A)  that  nevertheless  contains  all  the  essential  ingredi¬ 
ents  for  the  present  situation. 
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ANGLE  OF  INCIDENCE  (deg) 


NATURAL  REFLECTED  INTENSITIES 


Up  to  this  point,  we  have  been  considering  the  reflection  of  light  emitted  by  a  single  atom.  For  this 
case  the  cross  terms  in  (8),  those  proportional  to  the  product  of  Qs  and  Qp,  represent  the  effects  of 
interference  between  the  two  components  of  the  reflected  wave.  These  cross  terms  vanish  when  the 
incident  light  is  natural  for  the  following  reason.  A  natural  source  will  contain  many  atoms,  each 
radiating  independently  of  another  with  an  incident  amplitude  that  will  vary  about  some  average 
amplitude  and  an  angle  of  arrival  that  will  vary  at  random.  Furthermore,  the  amplitude  Eo  of  a  natural 
source  will  not  depend  on  the  angle  of  arrival,  a,  that  is,  a  natural  source  will  be  unpolarized.  If  we 
wish  to  apply  equations  (8)  to  a  natural  source,  then  we  must  average  (8)  over  all  angles  of  arrival. 
Since  the  averages  of  sin^  a  and  cos^  a  are  Vi  and  the  average  of  sin  2a  vanishes  (destroying  the 
interference  effects),  equations  (8)  for  the  reflection  of  natural  light  become 

_Nfy  R  2  Rp.  ^,2  /QN 

Rv  -  (cos/3) 

N^ky  Rs  Rp 

No  2^2 

where  the  reflected  radiance  A-'^is  assumed  to  be  proportional  to  E2*-E2  for  each  polarization  condi¬ 
tion  and  the  incident  sky  radiance  No  is  assumed  to  be  proportional  to  the  incident  intensity  Eo*'Eo. 

Equations  (9)  can  be  expected  to  hold  for  wave  numbers  up  to  and  including  those  in  the  infrared. 
However,  the  assumption  made  in  the  previous  paragraph  that  the  radiance  falling  on  the  sea  is  unpo- 
lanzed  will  fail  as  we  approach  the  visible  because  of  the  celebrated  effect  of  Rayleigh  scattering 
(McCartney,  1976),  dipole  scattering  from  irregularly  spaced  molecules  in  the  atmosphere  that,  in 
addition  to  introducing  a  blue  color  to  the  sky,  introduces  polarization  to  the  sky  radiance.  In  Appen¬ 
dix  D  we  estimate  that  the  total  polarization  of  the  sky  reaches  1  percent  at  a  wave  number  of  3500 
cm“*  (a  wavelength  of  2.9  ftm). 

It  is  evident  from  (9)  that 


Rf,+R,  =  R  (10) 

independent  of  the  value  of  p,  that  is,  independent  of  the  tilt  of  the  capillary  wave  facet.  This  implies 
that,  even  for  measurements  made  on  the  real  sea  containing  capillary  wave  facets  of  many  different 
slopes,  the  sum  of  the  measurements  with  the  polarizer  inserted  equals  the  measurement  with  the 
polarizer  removed  (assuming  a  perfect  polarizer). 

We  pause  here  to  emphasize  a  fact  obscured  by  our  notation,  namely  that  the  Fresnel  reflection 
coefficients,  and  Rp,  have  the  component  intensity  in  their  denominator  whereas  the  measured 
reflection  coefficients,  R^,  Ry,  and  R,  have  the  total  intensity  in  their  denominator. 
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NATURAL  EMITTED  INTENSITIES 


The  electromagnetic  theory  alluded  to  in  our  introductory  paragraph  shows  that  under  the  condi¬ 
tions  we  are  considering  in  this  problem,  plane  waves  incident  on  a  plane  boundary  between  two 
homogeneous  media,  the  equations  for  reflection  are  separable  into  two  independent  “5”  and  “p” 
components.  In  each  of  these  components,  the  reflection  coefficients  represent  the  ratio  of  reflected 
energy  to  incident  energy;  indeed  it  can  be  confirmed  by  inspection  of  the  complicated  equations  for 
Rs  and  Rp  (Stratton,  1941)  that  each  of  these  quantities  is  always  less  than  or  equal  to  one.  The 
remaining  energy,  the  energy  not  reflected  at  the  interface,  can  be  found  immediately  across  the 
interface  inside  the  other  medium  (water)  where  it  has  been  transmitted  by  the  boundary  and  where 
its  value  is  represented  by  Fresnel  transmission  coefficients  that  we  will  not  discuss.  The  transmitted 
energy  is  attenuated  by  absorption  as  it  propagates  away  from  the  interface. 

This  propagation  and  attenuation  might  continue  until  a  second  interface  was  encountered  where  a 
new  reflection  and  transmission  would  occur  very  similar  to  the  one  we  have  been  considering.  But 
now  we  are  involved  in  a  slightly  different  problem,  the  problem  of  reflection  by  a  thin  parallel  plate 
sample  of  water.  Let  us  return  to  our  original  case  for  which  the  water  is  thick;  none  of  the  energy 
can  escape  by  transmission;  and  all  the  energy  is  absorbed  within  the  water  itself.  Absorption  con¬ 
verts  the  optical  energy  into  heat;  the  water  is  heated;  the  black  body  emission  of  the  water  is 
increased;  and  we  are  implicitly  involved  in  another  new  development,  a  non-equilibrium  heating  of 
the  water  by  the  light.  This  process  would  continue  until  thermal  equilibrium  was  established  at  some 
higher  temperature  where  optical  absorption  was  balanced  by  thermal  emission.  This  balance  is 
expressed  by  Kirchoff’s  Law  (Bom  &  Wolf,  1970,  page  622),  which  says  that  for  each  wave  number 
(and,  we  assume  here,  for  each  component  “5”  and  “p”)  the  energy  emitted  by  the  water,  divided  by 
the  black  body  radiance  of  the  water  Nb^,  equals  the  fraction  A  that  the  water  absorbs  of  the  radiant 
energy  which  falls  upon  it: 


-^  =  A  =1. 


N, 


R, 


bb 


NT 

—^  =  a=\-r^ 

p  p 


N 


bb 


(11) 


The  Fresnel  emission  coefficients  As  and  Ap  are  shown  in  figures  4  and  5. 

Since  the  polarizer  acts  upon  these  emitted  intensities  in  the  same  way  that  it  acts  upon  reflected 
intensities,  the  horizontal,  vertical,  and  unpolarized  emitted  intensities.  Ah,  Ay,  and  A,  are  given  by  (9) 
with  Rs  replaced  by  A^,  and  Rp  replaced  by  A^: 


A  A 

A,=  ^  =  A(cosP)^+-^(sinP)^ 

^bh  ^  ^ 

=  A(sin  p)- +  ^(cos  p)^ 
A  A 

A= - =  A  +  = 

Nhb  2  2 


(12) 


10 


In  (12),  N^^%,  and  represent,  respectively,  the  emitted  radiance  detected  with  the  polarizer 
normal  to  the  plane  of  vision  (“horizontal”),  within  the  plane  of  vision  (“vertical”),  and  removed. 
Equations  (12)  show  that  the  horizontal  and  vertical  emissions  also  sum  to  the  unpolarized  emission 
in  the  real  sea  (for  a  perfect  polarizer). 

We  wish  to  comment  on  a  change  of  emphasis  brought  about  within  (12).  In  those  equations,  A 
and  R  are  energy  fractions  that  add  to  constant  values,  expressing  the  conservation  of  energy  (inci¬ 
dent  energy  is  absorbed  or  reflected).  However,  we  have  used  Kirchoff’s  Law  to  replace  absorbed 
energy  with  emitted  energy,  giving  A  and  R  different  denominators  (black  body  energy  and  incident 
energy,  respectively).  After  this  switch,  we  are  only  interested  in  the  numerators  by  themselves; 
hence,  there  is  nothing  to  be  gained  in  the  observation  that,  for  example,  Ay  +  =  Vi,  even  though 

this  happens  to  be  true. 
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INTEGRATION  OVER  SLOPES 


We  now  take  the  expressions  for  spectral  sea  radiance  we  have  derived  for  a  single-wave  slope, 
apply  a  weight  for  that  slope,  and  integrate  over  all  slopes.  Applying  this  to  (9)  and  (12)  gives: 

={NA}, 

Wf=(A'A>.  (13) 

N’‘>={N,Rl 


for  the  three  radiance  values  arising  from  sky  or  sun  reflections  and 


(14) 


for  the  three  radiance  values  arising  from  thermal  sea  emission. 

In  each  equation  of  (13)  and  (14),  the  angular  brackets  represent  the  following  integral  over  slope 
space  (Zeisse,  1995b): 


coar/2 
\Jr= const 


.  COSO)  ,,  ,,  ,,, 

•^cosG 


1/  cose 


co<7r/2 
U.  -const 


(15) 


where,  following  Cox  and  Munk  (1954),  p(^,  W)  stands  for  the  probability  that  a  wave  facet  will 
occur  with  a  slope  within  ±dl^/2  of  ^  and  ±dCy/l  of  when  the  wind  speed  is  W.  Cox  and  Munk 
(1954  &  1956)  obtained  an  expression  for  p,  the  occurrence  probability  density,  in  which  the  lowest 
order  term  is 


1 


-exp- 


2na 

a]  =0.000  + 3.16- 10'^  W 
gI  =0.003+ 1.92- 10’^  W 


1 

'~2 

(16) 


In  equation  (16),  a„  ^  and  Gc  ^  are  the  variances  in  ^  and  respectively,  and  W  is  the  wind  speed  in 
m  s“^. 
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The  factor  cos  (o/cos  0„  in  (15)  projects  the  horizontal  facet  area  toward  the  receiver,  weighting 
each  facet  by  the  amount  of  that  facet  “seen”  at  the  receiver.  This  weight  is  the  central  assumption  of 
a  sea-radiance  model,  contained  in  equations  (13)  through  (16),  which  has  been  more  fully  described 
elsewhere  (Zeisse,  1995c).  The  subscript  “r”  on  the  angular  bracket  serves  as  a  reminder  that  the 
receiver  direction  is  held  fixed  during  the  integration  over  sea  slopes  ^  and 

In  each  equation  of  (13),  No(Qs,  <|)5,  v)  represents  the  spectral  radiance  from  the  sky  or  sun  incident 
on  the  facet  from  source  coordinates  (6^,  that  are  related  to  the  (fixed)  receiver  coordinates  and 
the  ocean  slopes  by  the  law  of  specular  reflection.  When  the  radiance  is  from  the  sky,  the  integration 
is  over  all  slopes  in  the  sea  capable  of  reflecting  the  sky^  into  the  receiver.  When  the  radiance  is  from 
the  sun,  the  integration  is  only  over  the  tolerance  ellipse  (those  slopes  capable  of  reflecting  part  of 
the  solar  disk  toward  the  receiver).  In  (14),  Ntb  =  NbbiTsea,  v),  the  integration  is  over  all  slopes  in  the 
sea,  and  no  specular  reflection  is  considered  to  occur. 


^  Some  slopes  might  be  tilted  so  that  arriving  radiance  would  have  to  originate  below  the  horizon  (Gj  >  90°).  These  are 
excluded.  As  a  matter  of  fact,  shadowing  and  multiple  reflections  make  the  integration  unrealistic  when  radiance  arrives 
from  within  about  10°  of  the  horizon  (Gj  >  80°) 
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THE  SOLUTION 


The  previous  discussion  has  been  concentrated  solely  on  reflection  by  capillary  wave  facets  within 
the  footprint  of  a  single  pixel  from  an  image  of  the  ocean  surface.  Incident  values  have  been  those 
arriving  at  the  footprint,  and  reflected  values  have  been  those  immediately  leaving  the  footprint.  In 
most  practical  situations,  however,  the  receiver  is  located  a  significant  distance  away  from  the  foot¬ 
print,  typically  several  km.  Along  the  path  to  the  receiver,  the  reflected  and  emitted  radiance  are  atte¬ 
nuated,  and  path  radiance  is  accumulated. 

The  path  radiance  is  unpolarized  natural  light.  What  will  its  intensity  be  when  observed  through 
the  polarizer?  The  argument  follows  along  the  lines  of  those  employed  earlier.  If  we  consider  a  single 
atom  in  the  radiating  atmosphere,  then  we  may  describe  its  radiation  as  a  linearly  polarized  coherent 
plane  wave  propagating  along  the  optical  path  to  the  receiver.  The  incident  vector  amplitude,  E^,  will 
have  a  unique  direction  in  this  instance,  and  we  let  a  denote  the  angle  that  the  amplitude  makes 
with  the  axis  of  the  polarizer,  whatever  its  orientation.  The  polarizer  will  transmit  the  amplitude 
cos  a  corresponding  to  an  intensity  of  Eo*  Eo  cos^  a.  Averaging  over  all  values  of  a  to  represent  all 
atoms  radiating  in  this  natural  light  source,  we  conclude  that  a  polarizer,  if  perfect,  would  transmit 
half  of  the  unpolarized  value  regardless  of  its  orientation. 

Gathering  all  these  results  and  combining  reflected  and  emitted  intensities  in  each  direction,  we 
have  this  solution  for  the  sea  radiance  observed  with  the  polarizer  horizontal,  vertical,  and  removed: 


cos^  P  +  sin^  P ' 


+N^ 


N„  = 


sin^  P  +  cos^  P ' 


'l-/?^cos^  P-/?  sin^  p\  I  X 

bh\  )  r’'^/+ - — 

j^path  j, 


(17) 


N  =  N,+N^=^{N„-^—^ 


/  Rs+R. 


f  +  / 


Here  the  values  of  Rs  and  Rp  should  be  taken  from  equations  (A-9),  (A- 10),  and  (A-1 1),  and  equa¬ 
tion  (7)  should  be  used  for  p.  The  spectral  transmission  along  the  path  from  the  footprint  to  the 
receiver  is  represented  in  (17)  by  the  symbol  x.  The  function/has  been  included  in  (17)  to  represent 
the  relative  spectral  responsivity  of  the  instrument.  All  quantities  except  P  in  (17)  depend  on  the 
wave  number  v.  Broad-band  results  are  obtained  by  integrating  (17)  over  all  wave  numbers. 
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SEA  RADIANCE  RELATIONSHIPS 


It  is  convenient  to  separate  the  radiance  No  incident  on  each  facet  into  the  radiance  arriving  from 
the  sun  and  the  radiance  arriving  from  the  sky.  Then,  the  first  term  in  each  equation  of  (17)  can  be 
similarly  separated  into  one  part  due  to  the  sun  and  one  part  due  to  the  sky.  Denoting  these  by  the 
superscripts  ‘sun'"  and  ‘‘sky”  respectively,  we  can  express  (17)  in  more  compact  form  as 

={A^r + /  +{iv^"V2}/ 

A^v  =  { +  A/'v"''}x  /  +{ A^ ( 1 8) 

N 

Thus,  the  received  sea  radiance  consists  of  four  independent  components.  Let  us  refer  to  each  com¬ 
ponent  by  the  superscript;  where;  is  1  for  the  sun,  2  for  the  sky,  3  for  the  sea,  and  4  for  the  path. 
Total  radiance  will  be  referred  to  by  omitting  the  superscript  altogether.  We  define  the  degree  of 
polarization  of  the  ;th  component  to  be 


pj  =  Nl-Nj  _Ni-Nj 
N/l  +  NJ 


(19) 


Then,  given  the  total  radiance  and  polarization  of  each  component,  NJ  and  pj  respectively,  we  can 
find  all  other  quantities  from 


j 


p=^pj^ 

/  N 


(20) 


From  the  last  equation  we  see  that  the  total  polarization  consists  of  the  sum  of  the  polarizations  of 
each  component  weighted  by  the  fractional  radiance  contained  in  that  component. 
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DISCUSSION 


If  it  were  possible  to  observe  a  completely  flat  sea  at  a  zenith  angle  near  the  Brewster  angle,  about 
50  degrees,  reflected  light  would  be  completely  horizontally  polarized  because  of  the  vanishing  value 
of  Rp  in  the  plane  of  incidence,  which  is  vertical  in  this  case.  On  the  other  hand,  the  emitted  black 
body  radiation  would  be  completely  vertically  polarized  because  Ap  is  a  maximum  at  the  Brewster 
angle  (see  equations  (10)  and  figures  4  and  5).  If  reflection  were  to  dominate  emission,  or  vice  versa, 
then  a  high  degree  of  polarization  would  be  observed.  But  if  reflection  and  emission  were  compara¬ 
ble  in  intensity,  as  often  happens  in  the  long-wave  infrared  band,  each  would  tend  to  cancel  the  other 
giving  rise  to  rather  weak  polarization  of  the  radiance  from  a  perfectly  flat  sea. 

Such  an  ideal  condition  is  never  approached  in  the  real  sea,  however,  even  under  the  most  calm  of 
conditions.  There  is  always  a  breeze,  however  slight,  ruffling  the  surface  and  introducing  some  tilt  to 
the  capillary  wave  facets.  An  idea  of  how  much  tilt  there  is  can  be  gathered  from  (16),  which  states 
that  a  wind  speed  of  1  m  s“^  produces  a  distribution  of  wave  tilts  whose  standard  deviation  is  3 
degrees.  Although  this  is  small  compared  the  width  of  the  minimum  in  Rp  at  the  Brewster  angle,  it  is 
large  compared  to  the  angular  width  of  the  sun,  about  Va  of  a  degree.  Hence,  although  the  lightly 
ruffled  real  sea  under  calm  conditions  would  not  display  much  of  a  reduction  in  its  degree  of  polar¬ 
ization  below  that  of  a  totally  flat  sea,  the  amount  of  sun  glitter,  if  any,  would  be  greatly  reduced  by 
virtue  of  a  substantial  reduction  in  the  proportion  of  facets  within  the  tolerance  ellipse. 

As  the  wind  picks  up,  sky  reflection  and  thermal  emission  will  not  change  their  strength  substan¬ 
tially,  but  they  will  gradually  become  more  weakly  polarized  because  the  profusion  of  facet  tilts, 
involving  a  range  of  angles  of  incidence,  will  eliminate  the  sharp  distinction  between  “5”  and  “p” 
components  that  occurs  at  the  Brewster  angle.  The  wind  will  also  tilt  planes  of  incidence  as  it  tilts  the 
facets,  allowing  both  “j”  and  “p”  components  access  to  the  axis  of  “easy  passage”  though  the  polari¬ 
zer.  With  ever  higher  winds,  the  polarization  of  sky  reflection  and  sea  emission  will  gradually  be 
washed  away. 

On  the  other  hand,  the  polarization  of  the  glint  will  not  be  washed  away  as  the  wind  increases 
because,  regardless  of  wind  speed,  the  slopes  inside  the  ellipse  remain  exactly  the  same.  The  sun  is 
almost  a  point  source:  wave  slopes  very  close  to  the  specular  reflection  slope  (those  within  the  toler¬ 
ance  ellipse)  will  glint,  others  will  not.  From  another  point  of  view,  one  emphasized  recently  by 
Walker  (1994),  sun  glints  are  an  example  of  a  very  nonlinear  slope-to-radiance  transformation. 

Hence,  contrary  to  the  situation  with  sky  reflection  and  sea  emission,  the  polarization  of  sun  glint 
will  hardly  change  at  all  as  the  wind  increases. 

As  far  as  the  strength  of  sun  glint  is  concerned,  it  can  be  expected  to  go  through  a  maximum  as  the 
wind  increases.  The  reason  for  this  is  as  follows:  In  general,  the  tolerance  ellipse  will  be  located 
away  from  the  origin  of  slope  space  where  the  interaction  probability  density  is  small.  Hence,  the 
height  of  the  glint  column  (Zeisse,  1994a)  will  usually  be  small  under  calm  conditions.  As  the  wind 
picks  up,  the  interaction  probability  density  will  broaden,  increasing  the  height  of  the  glint  column. 
Finally,  at  high  winds,  the  flattening  of  the  interaction  probability  density,  which  accompanies  its 
broadening,  will  decrease  the  height  of  the  glint  column. 

Returning  to  the  subject  of  polarization,  path  effects  will  reduce  sea  polarization  for  two  reasons. 
First,  the  polarized  radiance  leaving  the  footprint  will  be  attenuated  by  the  path  transmission. 

Second,  unpolarized  radiance  will  accumulate  along  the  path  with  a  corresponding  dilution  in  the 
polarization  of  the  total  radiance. 
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In  figure  6,  the  radiance  of  each  component  is  shown  as  a  function  of  look  azimuth  with  respect  to 
the  azimuth  of  the  sun.  This  figure  is  taken  from  a  model^  of  the  real  sea.  The  spectral  region  is  the 
long-wave  band,  830  to  1250  cm"!  (8  to  12  pm)  at  low  spectral  resolution  (10  cm-^).  Navy  aerosols 
were  used  with  a  visibility  of  10  km  and  a  moderate  coastal  influence  (an  air  mass  parameter  ICSTL 
of  5  in  MODTRAN2).  The  solar  zenith  angle  is  70  degrees.  The  view  is  directly  into  a  wind  of 
10  ms-i  from  a  height  of  20  m  looking  down  onto  the  ocean  surface  at  an  angle  of  1  degree  with 
respect  to  the  horizon.  The  range  to  the  ocean  footprint  is  1 150  m  and  the  transmission  along  the 
path  to  the  footprint  is  82  percent.  These  are  windy  conditions  with  small  path  effects  and  strong  sun 
glints.  The  sky  polarization  is  +17.7  percent,  the  sea  polarization  is  -9.2  percent,  and  the  path  is 
assumed  to  be  unpolarized.  The  sun  polarization  varies  slightly  with  azimuth,  from  +22.7  percent  at 
an  azimuth  8°  away  from  the  sun  to  +27.6  percent  looking  directiy  along  the  sun’s  azimuth. 


Figure  6.  Sea  radiance  in  the  long-wave  band  for  a  rough  sea.  The  wind  speed 
is  10  m  s~\  the  solar  zenith  angle  is  70  degrees,  and  the  view  is  directly  up  wind 
from  a  height  of  20  m  looking  down  at  an  angle  of  1  degree  with  respect  to  the 
horizon.  Numbers  in  parentheses  refer  to  the  polarization  of  each  component. 


Figure  7  shows  the  radiance  in  the  mid-wave  spectral  band,  2000  to  3330  cm"!,  under  the  same 
sea  conditions  and  viewing  geometry.  The  transmission  along  the  path  to  the  footprint  in  this  case  is 
53  percent.  In  this  band,  the  path,  sea,  and  sky  radiance  lie  very  close  to  one  another;  for  clarity  only 
their  sum  has  been  plotted  in  figure  7.  The  combined  polarization  of  these  three  components  is  +5.1 
percent  using  equation  (19)  with  the  following  values  of  radiance  and  polarization  for  path,  sea, 
and  sky  respectively:  0.797  W  si^i  and  0  percent,  0.423  W  m-2  sr^  and  -13.9  percent,  and 


5  Figures  6, 7,  and  8  were  produced  by  a  polarized  version  of  “SeaRad,”  a  computer  code  for  the  prediction  of  sea 
radiance  (Zeisse,  1995c). 
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0.529  W  m  ^  si^*  and  28.2  percent.  The  polarization  of  the  sun  glint  varies  from  +28.9  percent  at  an 
azimuth  8°  away  from  the  sun  to  +35.3  percent  looking  directly  along  the  sun’s  azimuth. 


Figure  7.  Sea  radiance  in  the  mid-wave  band  for  a  rough  sea.  The  geometry 
is  the  same  as  for  the  previous  figure.  Numbers  in  parentheses  refer  to  the 
polarization  of  each  component. 


Figure  8  shows  the  total  polarization  in  each  band  for  these  conditions.  In  the  mid-wave  band,  sun 
glint  dominates  all  other  components,  and  the  total  radiance  along  the  center  of  the  pattern  is  prefer¬ 
entially  polarized  by  more  than  30  percent  in  the  horizontal  direction.  In  the  long-wave  band,  the 
solar  spectrum  is  weak  and  the  glints  have  diminished,  leaving  behind  competing  contributions  from 
sky  reflection,  sea  emission,  and  path  emission.  We  see  that  in  the  long-wave  band,  the  polarization 
of  the  total  radiance  is  negligible  for  the  case  chosen  here. 
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Figure  8.  Sea  polarization  in  the  mid-  and  long-wave  bands  for  a  rough  sea. 
The  conditions  are  the  same  as  for  figure  6.  A  negative  degree  of  polarization 
means  that  emission  dominates  reflection;  a  positive  degree  of  polarization 
means  that  reflection  dominates  emission.  Numbers  in  parentheses  refer  to 
the  equivalent  black  body  temperature  of  the  corresponding  signal. 


CONCLUSION 

For  unpolatized  incident  radiance,  horizontal  and  vertical  sea  radiance  sum  to  unpolarized  sea 
radiance. 

It  is  remarkable  that  this  relationship  holds  over  the  real  sea  since  it  reflects  and  emits  from  a 
bewildering  variety  of  wave  facet  tilts.  If  ‘The  sea  surface  is  a  mess”  (Tropf,  1990),  then  “The  polar¬ 
ized  sea  surface  is  a  bigger  mess,”  and  it  is  amazing  that  an  orderly  relation  such  as  this  one  has  man¬ 
aged  to  survive. 
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LIST  OF  SYMBOLS 


LATIN 


.4 

^s,p 

^h,v 

a 

b 

c 

f 

E 

k 

n 

No 

Nbb 

Nh,v 

N 

Ns.p 

Nh,v 

^s,p 

^h,v 

R 

U 

w 

X 

Y 

Z 


Unpolarized  emission  coefficient 
Fresnel  emission  coefficients 
Polarized  emission  coefficients 
Unit  vector  component  in  the  X  direction 
Unit  vector  component  in  the  Y  direction 
Unit  vector  component  in  the  Z  direction 
Relative  spectral  responsivity  of  receiver 
Amplitude  of  electromagnetic  vector 
Unit  vector  in  the  Z  direction. 

Complex  optical  index. 

Incident  spectral  radiance  (W  m~^  sr“*  (cm~*)“*) 

Spectral  black  body  radiance  of  the  sea 

Polarized  total  spectral  radiance 

Unpolarized  total  spectral  radiance 

Fresnel  radiance 

Polarized  radiance 

Fresnel  reflection  coefficients 

Polarized  reflection  coefficients 

Unpolarized  reflection  coefficient 

Unit  vector 

Wind  speed  (m  s~*) 

Coordinate  axis  in  up  wind  direction 
Coordinate  axis  in  cross  wind  direction 
Coordinate  axis  in  zenith  direction 
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GREEK 


a 

P 

5 

<}) 

V 

e 

p 

a 

X 

CO 

SUBSCRIPTS 

c 

h 

imag 

n 

o 

P 

r 

real 

s 

u 

V 
2 


Angle  of  arrival 

Angle  between  planes  of  incidence  and  vision 

Phase  of  plane  wave 

Azimuth  of  unit  vector 

Wave  number 

Zenith  angle  of  unit  vector 

Spectral  amplitude  reflectance  ratio 

Slope  variance 

Spectral  path  transmittance  (footprint  to  receiver) 
Angle  of  incidence 

Cross  wind 
Horizontal 

Imaginary  part  of  optical  index 
Facet  normal 
Incident  (sky  or  sun) 

Parallel  to  plane  of  incidence 
Receiver 

Real  part  of  optical  index 

Perpendicular  to  plane  of  incidence 

Up  wind 

Vertical 

Reflected 
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SUPERSCRIPTS 

path  Path  radiance  (footprint  to  receiver) 

sc  Scattered  solar  sky  radiance 

sea  Thermal  sea  radiance 

sky  Sky  radiance  reflected  by  the  sea 

sun  Solar  radiance  reflected  by  the  sea  (sun  glint) 

th  Thermal  sky  radiance 
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APPENDIX  A:  EQUATIONS  FOR  FRESNEL  REFLECTION  FROM  A 

CONDUCTING  SURFACE 

Let  electromagnetic  waves  of  frequency  v  propagate  through  a  medium  characterized  by  a  mag¬ 
netic  susceptibility  ji,  a  conductivity  o,  and  a  dielectric  constant  8.  Then  we  have  the  following  rela¬ 
tionships  (Stratton,  1941,  pp.  268  f¥): 


k  -a+ib 


03 

2n 

^=—^rea, 

C 

2% 

b  =  — n 

=  — n 

imag 

^0 

=  |ieG3  ^ 

2 


(A-1) 


where  is  called  the  separation  constant,  05  is  2%  times  the  frequency  of  the  wave,  Xq  is  the  free- 
space  wavelength,  and  rireai  and  ni^g  are  the  real  and  imaginary  parts,  respectively,  of  the  optical 
index  of  the  medium. 


Consider  the  case  in  which  there  is  Fresnel  reflection  at  the  boundary  between  two  media,  air 
(medium  2)  and  sea  water  (medium  1).  Assume  that  air  is  a  perfect  dielectric  (02  =  0)  and  that  sea 
water  is  a  conductor  (Oi  5^0).  Then  we  have  these  relationships  (Stratton,  1941,  page  501,  equation 
(48)): 


kl  =G5^e]|i,  +  /G5a,|i, 

A:j  =  a,  +  ffe, 


(A-2) 


for  the  water,  where  a\  and  b\  are  given  in  terms  of  the  material  constants  or  optical  indices  of  water 
by  (A-1),  and  these  relationships: 


^2  ®  ^21^2 

^2=^2 


(A-3) 


for  the  air. 

Let  CO  be  the  angle  of  incidence  and  define  p  and  q  (Stratton,  1941,  page  505,  equation  (69))  by 

=  al  -  b^  -  al  sin^  (O 

q^  +  -  -J -  al  sin^  (o)^ 

Then  (Stratton,  1941,  page  505,  equation  (74)) 


A-1 


p2_[^^2g-4i^2COSCof 

tan  8  = _ 2fij^2‘^2P^^^  ^ 

Ilf a^  COSTCO- 11 2 (q^ +  p^) 


and  (Stratton,  1941,  page  506,  equation  (77)) 


tan5„  = 


_  [M'2(3f-^f)cosQ)-p.,a29j  +[2|i2ajZ7,cosco-p,,fl2pf 
[\^2i^\  -  )  cos  (0+  |j,, 02^1  +[2|X2^J!A  cos  CO  +  |i,a2/^f 

2|X,p,2a2P^^  +p^-al  sin^  w)  cos  co 


^^«2  Kq^  +P')-P'l  W  +bf)  cos"  (0 


and  when  pi  =  P2  (Stratton,  1941,  page  506,  equation  (80)). 

_  Pp  _  [i? - «2  ®  co]^  +  p" 

pf  \q+a2  sin  CO  tan  co]^  +  p" 

tan5^tan(5  -S,)  = - 2a,PSino>tan  co _ 

^  of  sin"  CO  tan"  co  -  (9"  +  p" ) 


(A-5) 


(A-6) 


(A-7) 


When  a  and  b  are  expressed  in  terms  of  the  real  and  imaginary  parts  of  the  optical  index  of  the 
medium  using  (A-1),  the  factor  2nAo  cancels  from  each  of  equations  (A-5)  through  (A-7).  Hence 
in  these  equations,  23i/Xo  can  be  set  equal  to  one,  and  the  real  part  of  the  index  can  be  used  in  place 
of  a,  and  the  imaginary  part  of  the  index  can  be  used  in  place  of  b.  With  these  substitutions,  and  with 


^Ireal  ~  ^real 
\imag  imag 


(A-8) 


the  equations  for  q  and  p  become 


q  -P  =n^eal  -  n,nu,g  "  Sm  CO 


q^  +  P^  =  ^  +  («L  -  -  sin"  CO^ 


the  perpendicular  reflection  coefficient  and  phase  become 

2_[g-cos  cof +  /  ^ 

[^  +  costof-fp"  ^ 

~  2pcosto 

tano,  = - 7 - . 

cos  to  -  (<?"  +  p") 


(A-9) 


(A-10) 


A-2 


the  parallel  reflection  coefficient  and  phase  become 


and 


[fc/  -«L,)cos  (P-q 

1 

2«..«A„„,cos0-p]' 

[(«L  -"L,)cos  <P+q 

2 

+ 

2«„aA™,cos0-hp]' 

2p(g^  +  p^-  sin^  0)  cos  0 
tano  „  = - - 

q^  +  p^-  «,,+n,^,Jcos'0 


(A-11) 


P 


2 


Pp  _  [^-sin0  tan0]^  +  p^ 
Ps  +  sin0  tan0]^ 


tan6^tan(5  -5)  = _ , 

’’  '  sin^ 0  tan^ (O  -{q^  +  p^) 


(A-12) 


When  one  considers  that  the  quantities  q  and  p  are  functions  of  the  angle  of  incidence  as  well  as 
of  all  the  parameters  of  the  medium,  the  complexity  of  what  appeared  at  first  to  be  the  simplest  of 
problems — the  reflection  of  a  plane  wave  from  a  plane,  absorbing  surface — is  truly  amazing.  The 
formulas  are  far  too  involved  to  be  understood  from  a  casual  inspection,  and  the  nature  of  the  reflec¬ 
tion  phenomena  becomes  apparent  only  when  one  treats  certain  limiting  cases.”  (Stratton,  1941,  page 
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APPENDIX  B:  FORTRAN  CALCULATION  OF  FRESNEL  REFLECTION 

FROM  SEA  WATER 


0  SUBROUTINE  Rho( omega,  v,  Rs,  Rp) 

CU  USES  SeaData 

COMMON  /Sealndex/  alphaOl (100 ) ,  alpha02(20), 
+  betaOl  (100),  beta02  (20) 


*  Calculates  reflectivity  of  sea  water  between  52.63  cm-1  and  * 

*  25,000  cm-1  using  equations  (74)  and  (78)  of  Stratton,  “Electro-  * 

*  magnetic  Theory,"  1941,  page  505,  ff.  The  sea  water  is  assumed  to  * 

*  be  a  conducting  mediiim;  both  real  and  imaginary  parts  of  the  * 

*  index  of  water  are  used.  The  notation  of  Stratton  is  adhered  to  * 

*  as  closely  as  possible  except  that  I  have  set  alpha  equal  to  the  * 


*  real  part  of  the  index  [rather  than  (2*pi) /lambda*Real (n) ]  and  beta  * 


*  equal  to  the  imaginary  part  of  the  index,  [rather  than  * 

*  (2*pi) /lambda*Imag(n) ] ,  which  can  be  done  because  the  factor  * 

*  (2*pi) /lambda  cancels  out  in  these  equations.  Lambda  is  the  free  * 

*  Space  wavelength.  ★ 

* 


*  The  calculation  is  for  a  single  flat  wave  facet,  sea  water  on  one  * 

*  side  and  air  on  the  other  side.  * 


*  Omega  is  the  angle  of  incidence  in  radians,  and  v  is  the  wave  * 

*  number  in  cm-1;  these  are  both  inputs.  The  outputs  are  Rp  and  Rs .  * 

*  Rp  is  the  reflectivity  that  would  be  observed  with  a  polarizer  * 

*  oriented  parallel  to  the  plane  of  incidence.  Rs  is  the  reflectivity  * 

*  that  would  be  observed  with  a  polarizer  oriented  perpendicular  to  * 

*  the  plane  of  incidence.  To  find  the  unpolarized  reflectivity,  take  * 

*  (Rp  +  Rs)/2. 

* 

★ 

*  Last  revision:  December  28,  1995  * 
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The  four-point  interpolation  functions  are: 

WM(x)  =  (x  -  l.)*(x  -  2.)*x/6. 

W0(x)  =  (x  -  l.)*(x  -  2.)*(x  +  l.)/2. 

Wl(x)  =  (x  +  l.)*(x  -  2.)*x/2. 

W2 (x)  =  (x  +  l.)*(x  -  l.)*x/6. 

If  ((omega  .It.  0.)  .or.  (omega  .gt.  1.57080))  then 
omega  is  out  of  bounds;  set  reflectivities  to  0%  and  return 
Rp  =  0. 

Rs  =  0. 

Return 
End  If 

If  (v  .eq.  0.)  then 

set  reflectivities  to  100%  and  return: 

Rp  =  1. 

Rs  =  1. 

Return 
End  If 
W  =  1.E4/V 

If  (W  .LT.  0.399999)  then 
print  error  message  and  return: 

Rp  =  0. 

Rs  =  0 . 

Write  (6,  1000)  v 
Return 
End  If 

If  (0.4  .le.  W  .and.  W  .le,  19.8)  then 

use  Lagrange  4  point  interpolation  on  Block  01  data  which 
are  at  0 . 2  um  spacing  between  0.2  and  20  um: 

I  =  Int(W/0.2) 

Fr  =  Mod(W,  0.2) /O. 2 

alphal  =  W2 (Fr) *alpha01 (I  +  2)  -  W1 (Fr ) *alpha01 ( I  +  1) 

+  WO (Fr) *alpha01 (I)  -  WM (Fr ) *alpha01 ( I  -  1) 

betal  =  W2 (Fr) *beta01  (1+2)  -  W1 (Fr ) *beta01  (I  +  1) 


+  wo (Fr) *beta01  (I) 


-  WM(Fr) *beta01  (I  -  1) 


End  If 

If  (19.8  .It.  W  .and.  W  .It.  190.)  then 
C  use  Lagrange  4  point  interpolation  on  Block  02  data  which 

C  are  at  10  um  spacing  between  10  and  200  urn: 

I  =  Int(W/10.) 

Fr  =  Mod(W,  10. ) /lO. 

alphal  =  W2 (Fr) *alpha02 (I  +  2)  -  W1  (Fr ) *alpha02 (I  +  1) 

+  +  WO  (Fr) *alpha02 (I)  -  WM (Fr ) *alpha02 ( I  ~  1) 

betal  =  W2 (Fr) *beta02  (1+2)  -  W1 (Fr ) *beta02  (I  +  1) 

+  +  WO (Fr) *beta02  (I)  -  WM (Fr ) *beta02  (I  -  1) 

End  If 

If  (190.  .LE.  W)  then 
C  print  error  message  and  return: 

Rp  =  0. 

Rs  =  0. 

Write  (6,  1000)  v 
Return 
End  If 

g  =  alphal* *2  ~  betal* *2  -  sin (omega) * *2 
h  =  4 . * {alphal**2 ) * (betal**2 ) 
p  =  sqrt(0.5*(~g  +  sqrt(h  +  g**2))) 
q  =  sqrt(0.5*(+g  +  sqrt(h  +  g**2))) 

C  Stratton,  Equation  (74),  p.  505: 

c  =  (q  -  cos (omega) )* *2  +  p**2 

d  =  (q  +  cos (omega) )* *2  +  p**2 

Rs  =  c/d 

C  Stratton,  Equation  (77),  p.  506: 

e  =  (  (alphal* *2  -  betal * *2 ) *cos (omega)  -  q)**2 

f  =  ({alphal**2  -  betal**2 ) *cos (omega)  +  q)**2 

t  =  (2*alphal*betal*cos (omega)  -  p)**2 

u  =  (2*alphal*betal*cos (omega)  +  p)**2 

Rp=  (e+t)/(f+u) 

Return 
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1000  FORMAT  ('  *****  WARNING  -  INPUT  FREQUENCY  =  1PG12.5, 

+  /,  '  OUTSIDE  VALID  RANGE  OF  52.63  TO  25,000  CM-1  ******> ^/) 

End 

BLOCK  DATA  SeaData 

COMMON  /Sealndex/  alphaOl ( 100 ) ,  alpha02(20), 

+  betaOl  (100),  beta02  (20) 


*  These  data  for  the  optical  index  of  water  have  been  taken  from  * 

^  the  literature.  From  0 . 2  to  2  microns  (blocks  01  up  to  second  * 

*  entry  of  row  B)  and  from  10  to  200  microns  (blocks  02)  the  data  * 

*  are  from  G.  M.  Hale  and  M.  R.  Querry,  "Optical  Constants  of  Water  * 

*  in  the  200-nm  to  200-um  Wavelength  Region,"  Appl .  Opt.  3,  555  * 

*  (1973) .  These  data  are  for  pure  water.  * 

*  * 

*  From  2.2  to  20  microns  (blocks  01  from  the  third  entry  of  row  B  * 

*  to  the  end),  the  data  are  from  M.  R.  Querry,  W.  E.  Holland,  R.  C.  * 

*  Waring,  L.  M.  Earls,  and  M.  D.  Querry,  "Relative  Reflectance  and  * 

*  Complex  Refractive  Index  in  the  Infrared  for  Saline  Environmental 

*  Waters,"  J.  Geophys .  Res.  82,  p.  1425  (1977),  Table  3,  Pacific  * 

*  Ocean  columns.  These  data  are  for  salt  water.  * 


C  Real  part  of  the  index  of  sea  water  from  0.2  to  20  microns  in 

C  steps  of  0.2  microns: 


DATA 

alphaOl 

/ 

A 

1.396, 

1.339, 

1.332, 

1.329, 

1.327, 

1.324, 

1.321, 

1.317, 

B 

1.312, 

1.306, 

1.303, 

1.287, 

1.251, 

1.151, 

1.384, 

1.479, 

C 

1.421, 

1.388, 

1.368, 

1.355, 

1.347, 

1.339, 

1.335, 

1.335, 

D 

1.332, 

1.324, 

1.312, 

1.296, 

1.268, 

1.271, 

1.371, 

1.353, 

E 

1.340, 

1.330, 

1.324, 

1.319, 

1.314, 

1.307, 

1.302, 

1.297, 

F 

1.291, 

1.286, 

1.279, 

1.272, 

1.268, 

1.264, 

1.258, 

1.249, 

G 

1.240, 

1.229, 

1.218, 

1.204, 

1.190, 

1.177, 

1.165, 

1.151, 

H 

1.140, 

1.132, 

1.124, 

1.119, 

1.121, 

1.126, 

1.134, 

1.142, 
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I 

1.152, 

1.164,  1.177,  1.189, 

1.201, 

1.212, 

1.224, 

1.234 

J 

1.242, 

1.253,  1.261,  1.273, 

1.284, 

1.296, 

1.309, 

1.320 

K 

1.331, 

1.339,  1.349,  1.358, 

1.366, 

1.379, 

1.390, 

1.399 

L 

1.408, 

1.417,  1.426,  1.435, 

1.443, 

1.450, 

1.458, 

1.464 

M 

1.470, 

1.474,  1.477,  1.480 

/ 

C  R03I  p03rt  of  til©  in<i©x  of  s©fi  wstGX"  firoiti  10  to  200  iniciroris  in 
C  steps  of  10  microns: 

DATA  alpha 02  / 

A  1.218,  1.480,  1.551,  1.519,  1.587,  1.703,  1.821,  1.886, 

B  1.924,  1.957,  1.966,  2.004,  2.036,  2.056,  2.069,  2.081, 

C  2.094,  2.107,  2.119,  2.130  / 

C  Imaginary  part  of  the  index  of  sea  water  from  0.2  to  20  microns  in 

C  steps  of  0.2  microns : 

DATA  betaOl  / 


A 

0 

.000, 

0 

.000, 

0 

.000, 

0 

.000, 

0 

.000, 

0 

.000, 

0 

.000, 

0 

.000, 

B 

0 

.000, 

0 

.001, 

0 

.000, 

0 

.001, 

0 

.003, 

0 

.114, 

0 

.263, 

0 

.085, 

C 

0 

.018, 

0 

.005, 

0 

.003, 

0 

.004, 

0 

.007, 

0 

.011, 

0 

-016, 

0. 

.016, 

D 

0. 

.013, 

0 

.011, 

0 

.010, 

0 

.013, 

0 

.032, 

0 

.108, 

0 

.087, 

0, 

.044, 

E 

0, 

.035, 

0, 

.033, 

0 

.032, 

0, 

.031, 

0, 

.031, 

0, 

.032, 

0 

.032, 

0. 

.033, 

F 

0. 

.034, 

0, 

.036, 

0 

.038, 

0. 

.041, 

0, 

.044, 

0. 

.046, 

0, 

.046, 

0. 

.047, 

G 

0. 

.048, 

0. 

.050, 

0. 

.054, 

0. 

.060, 

0. 

.068, 

0. 

.079, 

0. 

.091, 

0. 

.107, 

H 

0. 

.125, 

0. 

.145, 

0. 

.166, 

0. 

.191, 

0. 

.216, 

0. 

.239, 

0. 

.260, 

0. 

.279, 

I 

0. 

.297, 

0. 

.313, 

0. 

.326, 

0. 

.338, 

0. 

.347, 

0, 

.357, 

0. 

.363, 

0. 

.371, 

J 

0. 

,377, 

0, 

.385, 

0. 

.393, 

0. 

.401, 

0. 

.407, 

0. 

.413, 

0. 

.417, 

0. 

.418, 

K 

0. 

,420, 

0. 

,422, 

0. 

.424, 

0, 

.427, 

0, 

.430, 

0. 

.432, 

0. 

,432, 

0. 

,432, 

L 

0. 

431, 

0. 

430, 

0, 

.429, 

0. 

.427, 

0. 

,425, 

0. 

,423, 

0. 

,420, 

0. 

416, 

M 

0. 

414, 

0. 

410, 

0. 

.406, 

0. 

393 

/ 

C  Imaginary  part  of  the  index  of  sea  water  from  10  to  200  microns  in 

C  steps  of  10  microns: 

DATA  beta02  / 

A  0.051,  0.393,  0.328,  0.385,  0.514,  0.587,  0.576,  0.547, 

B  0.536,  0.532,  0.531,  0.526,  0.514,  0.500,  0.495,  0.496, 

C  0.497,  0.499,  0.501,  0.504  / 

End 
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APPENDIX  C:  THE  ANGLE  BETWEEN  THE  PLANES  OF  INCIDENCE 

AND  VISION 


In  this  appendix  we  derive  an  expression  for  p,  the  angle  between  the  plane  of  incidence  and  the 
plane  of  vision.  Let  Up  be  a  unit  vector  normal  to  U;.  pointing  along  the  direction  of  E^p  in  the  plane 
of  incidence  and  let  Uv  be  a  unit  vector  normal  to  U;-  pointing  along  the  vertical  axis  of  polarization 
in  the  plane  of  vision.  Then  (see  figure  3) 

cosp=U^-U,  (C-1) 

To  find  P  we  must  find  the  coordinates  of  Up  and  Uv  and  then  take  their  inner  product. 

First  we  will  find  the  coordinates  of  Up.  Since  it  is  normal  to  U^- 

U,-U,=0  (C-2) 

and  since  it  is  in  the  plane  of  incidence,  it  can  be  written  as  a  linear  combination  of  two  other  vectors 
in  that  plane: 


U  =AU  +5U 


(C-3) 


Setting  the  inner  product  of  (C-3)  with  U;-  equal  to  zero,  recalling  that  co  is  the  angle  between  U„ 
and  Ur,  and  normalizing  the  length  of  Up  to  one,  we  find  that 


±U  = 


U„  -  COSO) 
sin  CO 


(C-4) 


Next  we  want  the  coordinates  of  Uv  These  might  appear  to  be  trivially  available  since  Uv  is  “verti¬ 
cal”  and  coincident  with  the  unit  vector  k  along  the  Z  axis  with  coordinates  (0,  0,  1).  However,  Uv  is 
not  vertical  since  it  is  normal  to  U^,  a  vector  that  is  generally  not  horizontal.  The  coordinates  of  Uv 
must  be  found  by  considering  how  the  polarizer  is  actually  adjusted  in  practice:  first,  the  polarizer  is 
set  horizontal  (no  Z  component)  and,  second,  the  polarizer  is  turned  by  90  degrees.  So,  to  find  the 
coordinates  of  Uv,  we  have  to  find  the  coordinates  of  U/,  and  then  construct  Uv  normal  to  U/^  and  U^. 
Since  U/,  is  normal  to  U;- 


u.-  U  =0 


(C-5) 


and  since  U/,  is  horizontal 


U,-k  =  0.  (C-6) 

These  conditions,  plus  the  requirement  that  the  length  of  U/,  be  one  result,  after  some  algebra,  in  the 
coordinates 


±U;^  =(sin(j)^,  -cos<j)^,  0)  (C-7) 

The  calculation  for  Uy  follows  similar  lines  but  will  be  avoided  altogether  by  noting  that,  as  shown 
in  figure  3,  P  is  the  complement  of  the  angle  between  Up  and  Vh : 


C-1 


(C-8) 


sinp  =  cos  (7t/2  -  p)  =  U  •  ^ 

sin  CO 

where  we  have  used  the  fact  that  U/,  is  normal  to  U^.  According  to  (1)  the  coordinates  of  U„ 

U„  =(sin0„cos(^„,  sine„sin<|)„,  cosGj 

and  substitution  of  (C-9)  and  (C-7)  into  (C-8)  gives 

.  -  sin0„sin((j)  -<j)  ) 

sin  p  = - 2 - 

sin  CO 

for  the  angular  separation  between  the  planes  of  incidence  and  vision. 


are 

(C-9) 


(C-10) 
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APPENDIX  D:  POLARIZATION  OF  THE  INFRARED  SKY 


The  purpose  of  this  appendix  is  to  estimate  the  wave  numbers  at  which  the  total  polarization  of  the 
infrared  sky  reaches  1  percent  and  10  percent.  The  reader  is  referred  to  (McCartney,  1976)  for  a 
more  complete  discussion  of  scattered  solar  radiance. 


In  the  infrared  the  primary  source  of  sky  radiance  is  thermal  emission  of  the  atmosphere.  In  the 
visible  the  primary  source  of  sky  radiance  is  scattered  sunlight.  We  will  assume  that  for  all  wave 
numbers  v  throughout  the  infrared  and  visible,  the  total  spectral  sky  radiance  No(\)  consists  of  a 
combination  of  spectral  thermal  radiance  N%(v)  and  spectral  scattered  solar  radiance  N^%(v).  Then 
equation  (20)  gives  for  the  total  sky  polarization 


p  ^  C(v)A^r(v)+r(v)A^f(v)  , 

Ar-(v)+iVf(v)  ^ 

where  P^%(v)  is  the  polarization  of  the  scattered  solar  radiance  and  P%(v)  is  the  polarization  of  the 
thermal  radiance.  The  polarization  of  the  scattered  solar  radiance  varies  between  —10  percent  and 
+90  percent  depending  on  the  observation  direction  (McCartney,  1976,  page  212),  and  for  the  pur¬ 
poses  of  this  estimate  we  shall  assume  that  it  is  always  +100  percent.  The  polarization  of  the  thermal 
radiance  vanishes  since  thermal  radiance  is  a  natural  source.  Hence  the  problem  posed  above 
becomes  the  problem  of  finding  the  values  of  v  for  which 


iVf(v) 

A^“(v)  +  N'J'iv) 


(D-2) 


is  equal  to  1  percent  and  10  percent. 

The  thermal  spectral  radiance  can  be  approximated  by  the  spectral  radiance  of  a  black  body  at 
absolute  temperature  300  K: 


N'J"  ^  (D-3) 

The  following  expression  has  been  given  (McCartney,  1976,  p.  208)  for  single  scattered  solar 
radiance  found  at  the  bottom  of  a  plane  parallel  atmosphere  when  the  sunlight  arriving  at  the  top  of 
the  atmosphere  has  a  spectral  irradiance  Eo(v): 


N^f(v ,  e) = 

Pm 

In  this  expression,  9  is  the  supplement  of  the  angle  between  the  direction  to  the  sun  and  the  direction 
to  the  point  of  observation  in  the  sky;  is  the  zenith  angle  of  the  solar  direction;  Q,  is  the  zenith 
angle  of  the  observation  direction;  and  Zg  is  the  reduced  height  of  the  atmosphere.  The  volume 
angular  scattering  coefficient,  ^^(9),  is  the  intensity  scattered  into  an  infinitesimal  solid  angle  da>  at 
an  observation  angle  9  with  respect  to  the  direction  of  the  incident  light.  For  unpolarized  incident 
light,  it  is  given  by 


exp  (-  P  sec  C  ^)  -  exp  (-  P  sec  ^  ^  ) 

l-secC,cosC„ 


(D-4) 
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P.(6)  = 


(l+cos^S) 


(D-5) 


where  n  is  the  refractive  index  of  the  molecular  atmosphere  consisting  of  dipoles  per  unit  volume. 
Finally,  the  total  volume  scattering  coefficient 


rm  Jq 

is  the  ratio  of  the  flux  scattered  in  all  directions,  by  unit  volume  of  gas,  to  the  irradiance  of  the  inci¬ 
dent  flux'  .  The  use  of  a  plane  parallel  atmosphere  causes  errors  of  less  than  2  percent  for  all  zenith 
angles  ^  less  than  80  degrees  (sec  ^  less  than  6). 

We  wish  to  show  that  the  argument  of  each  exponential  in  (D-4)  is  small.  Values  of  (n^  -  1)2  are 
weakly  dependent  on  wave  number  and  stay  within  several  percent  of  3  x  10-2  throughout  the 
infrared  and  visible  (McCartney  1976,  page  348).  N  is  given  by  the  number  of  molecules  per  unit 
volume  in  the  uniform  reduced  atmosphere  at  standard  temperature  and  pressure,  namely,  2.547  x 
10'9  cm-3.  This  means  that  the  value  of  %n  is  approximately  10-2  km-'  at  a  wave  number  of  10,000 
cm-'.  Next,  the  reduced  height  of  the  1962  Standard  Atmosphere  is  8.00  km  (McCartney,  1976,  page 
92),  and  sec  ^  is  less  than  6,  so  the  arguments  of  the  two  exponentials  in  (D-4)  are  indeed  less  than 
one  in  the  infrared  and  we  can  approximate  (D-4)  by 


sec  C,- sec  C, 

1  -  sec  ^  ^  cos  ^  p 

Furthermore,  we  can  replace  (1  +  cos2  9)  in  (D-5)  by  its  maximum  value  of  2,  neglect  the  square 
bracket  in  (D-7)  whose  value  is  on  the  order  of  unity,  and  express  the  solar  irradiance  at  the  top  of 
the  atmosphere  in  terms  of  the  solar  radiance  there  times  the  solid  angle  cq,  subtended  by  the  sun, 
which  is  6.763  x  10-^  sr.  We  thus  arrive  at 


Af“(v,e)=£ri3.(e)z: 


JVr(v)“JV«(5900i:,v)-a). 


6.29  •  10-^^  (5900  K,\  )  -v 


for  the  spectral  radiance  of  a  plane  parallel  Rayleigh  sky  at  the  surface  of  the  earth  in  the  infrared 
visible  spectral  regions.  In  (D-8)  we  have  gone  on  to  assume  that  the  solar  spectral  radiance  is  given 
by  Planck’s  equation  for  a  black  body  NibiZ  v)  at  absolute  temperature  T=  5900  K. 


*  In  expressions  (D-5)  and  (D-6),  it  is  helpfiil  to  remember  that  n^  -  1  is  closely  proportional  to  N  so  that  the  scattering 
coefficients  are  proportional  to  the  number  of  scattering  centers  N  per  unit  volume  as  expected. 
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Table  D-1  lists  some  values  of  spectral  sky  radiance  according  to  (D-8)  and  its  polarization  accord¬ 
ing  to  (D-2); 


Table  D-1 .  Spectral  sky  radiance  and  polarization. 


V 

Nbb(5900  K,  V) 

No®‘=(v) 

Nbb(300  K,  v) 

P(v) 

(cm*'') 

(W  m-2  sri  (cm--*)-!) 

(W  m~2  sr^  (cm“i)~'') 

(W  m-2  sH  (cm-'')-^) 

(%) 

3,000 

269 

1. 37  X  10-7 

1.64  X  10-4 

0.1 

3,500 

341 

3.22  X  10-7 

2.37  X  10-5 

1.3 

4,000 

415 

6.68  X  10-7 

3.23  X  10-6 

17.1 

4,500 

489 

1.26  X  10-6 

4.18  X  10-7 

75.1 

5,000 

562 

2.21  xlO-6 

5.22  X  10-® 

97.7 

Hence  the  required  estimate  is 


P(v)«l%arv  «  3,500  cm-'  (2.9  |Llm) 
P(v)«10%  arv  «  4,000  cm-'  (2.5  pm) 
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